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Abstract 

' We analyse the propagation of electromagnetic waves in magnetoelectric media. Recently, Feigel 

has predicted that such a medium may "extract momentum from vacuum" in the sense that 

o : 

the total momentum of the virtual waves (vacuum fluctuations of the electromagnetic field) is 

O 

nontrivial. Our aim is to check the feasibility of this effect. The crucial point in our study is 

>V 

an assumption of the finite size of the magnetoelectric sample, which allows us to reduce the 

O-i. 

calculation of the momenta and forces of the electromagnetic waves acting on the sample to the 
vacuum region outside of the medium. In this framework, we demonstrate that, in contrast to 
Feigel, the total force caused by the virtual is zero, with an appropriate count of the modes that 
■ should be taken into account in this effect. Furthermore, we find that the two irreducible parts 

of the magnetoelectric matrix behave differently in the possible Feigel effect. Going beyond the 
original scheme of the virtual electromagnetic waves, we propose an experimental scheme which is 
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suitable for the measurement of the magnetoelectric susceptibilities of the medium with the help 
of real electromagnetic waves. 
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I. INTRODUCTION 



Phenomenological macroscopic electrodynamics is based on the well known experimental 
observations that an external electric field can induce a polarization of a medium, whereas 
an external magnetic field can induce a magnetization of matter. As a result, the electric 
and magnetic excitations (T>, 7i) (comprising a 2-form H) are functions of the electric and 
magnetic field strengths (E, B) (collected in a 2-form F) and of the permittivity e and the 
permeability of the medium. In the simplest case of an isotropic medium at rest, the 
constitutive relations read V = ee E and H = (/i/i ) _1 -B. Here e and fi are the electric 
and magnetic constants (permittivity and permeability of the vacuum). However, in the 
1960s it was theoretically predicted [1] and experimentally confirmed j^, 3] that certain 
media become electrically polarized when placed into a magnetic field or are magnetized 
when put into an electric field. Since then, such a magnetoelectric effect was observed 
for many substances and studied in great detail both theoretically and experimentally, see 
for the reviews {4J, 5], for example. The constitutive relation of a medium at rest is then 
modified to V — ee E + (3 ■ B and 7i = (/i/i ) _1 5 — (3 T ■ E, where the traceless 3x3 
matrix {3 describes the magnetoelectric properties ( T denotes the transposed matrix). The 
magnetoelectric effect can be observed in a certain class of media, but it also can be induced 
in media which are put into external electric and magnetic fields. 

Recently, attention to magnetoelectric media was attracted in connection with an inter- 
esting new effect predicted by Feigel He noticed that the propagation of electromagnetic 
waves in a magnetoelectric medium is essentially asymmetric in the sense that the waves 
moving in opposite directions carry different momenta. Then, calculating the total momen- 
tum of the virtual waves (or "vacuum fluctuations" of the electromagnetic field) which are 
present in the sample, he concluded that this quantity is nontrivial. In this sense, the Feigel 
effect predicts the extraction of momentum from vacuum. This possibility was discussed 
in [3, [a, [9, [h], H, 12)]. In particular, in [3| a certain similarity of the Feigel effect with 
the Casimir effect was noticed, and the computation of the total momentum transfered by 
the virtual waves was attacked by means of the Green's function method. In this paper, 
we analyse the Feigel effect in a sample of finite size which was somehow neglected in the 
previous works. 

It seems clear that the discussion of the Feigel effect is related to the definition of the 
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energy and momentum of the electromagnetic field in a medium,and moreover, in a mov- 
ing medium (since a nontrivial velocity of a sample is predicted [6|). However, as it is well 
known, the issue of the energy and the momentum of the electromagnetic field (and of waves, 
in particular) in dielectric and magnetic media has a long and controversial history. The 
discussion of this issue began with the investigations of Min 
instein and Laub [15(|. The problem is reviewed in 16 
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13), Abraham 14|, and 
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em was not settled neither theoretica 



23) . e.g.). In the recent paper |24( we have proposed a 



ly, nor experimentally 



consistent definition of the electromagnetic energy-momentum in media, and furthermore, a 
variational approach was developed in [25j| for the moving magnetoelectric medium specifi- 
cally for the study of the Feigel effect. An interesting technical observation is then that the 
magnetoelectric matrix (3 induces a term in the energy-momentum tensor which describes 
an additional flux of energy and momentum. Such an additional term is present even in the 
medium at rest, and it vanishes when the magnetoelectric properties are absent. This fact 
lends some support to the feasibility of the Feigel effect, at least on a qualitative level. 

Here, however, we will analyse the possibility of the Feigel effect avoiding the problem of 
the electromagnetic energy-momentum in the medium. The crucial point is the assumption 
of a finite size of the magnetoelectric sample. Then we notice that the virtual waves, 
which are excited in the medium, are not confined to the sample. Since the boundaries 
are not impenetrable, the vacuum fluctuations exist everywhere, inside the medium as well 
as outside of it. Since the virtual electromagnetic waves in these regions of space are related 
by the jump conditions across the boundaries, we can eventually replace the evaluation of 
the electromagnetic momentum and forces inside the medium by the computation of these 
quantities in free space. In vacuum the energy-momentum is uniquely defined, and this 
simplifies the analysis of the possible Feigel effect to a considerable extent. 

The structure of the paper is as follows. Sec. HTlpresents some general material, introduces 
the basic notions and fixes the notation. In Sec. II I II we analyse the constitutive relation of 
the magnetoelectric medium. In particular, we provide an irreducible decomposition of 
the magnetoelectric matrix, which proves to be convenient for the subsequent theoretical 
analysis. Sec. [IV] demonstrates a general feature of the wave propagation in magnetoelectric 
media, namely the birefringence. We demonstrate that for a medium characterized by the 
second irreducible piece, the generic Fresnel covector surface factorizes into a product of 
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two light cones. We find the two corresponding optical metrics. In Sec. |V]the propagation 
of plane electromagnetic waves is studied in a magnetoelectric medium of the Feigel type. 
From the jump conditions on the boundaries between the vacuum regions and the medium, 
we derive relations between the transmitted and the reflected waves. These relations are 
then used in Sec. I VI I for the calculation of the electromagnetic energy, the momentum and 
the force acting on the magnetoelectric sample. 

II. SOME BACKGROUND MATERIAL 



If the 4-dimensional (4D) electromagnetic excitation 2-form is denoted by 36| 



H =- H A da + V = -Hijdx* A dx j (2.1) 
and the 4D electromagnetic field strength 2-form by 

F = E Ada + B = ^F^dx 1 A dx j , (2.2) 
then the Maxwell equations read 

dH = J, dF = 0. (2.3) 

Here a is a time parameter and J=—j/\da + p the 4D electric current. The 3D magnetic 
excitation Ti. = 7i a dx a is a 1-form, the electric excitation T> = \T> a bdx a A dx b a 2-form. 
Analogously, we have the 3D electric field strength E = E a dx a and the 3D magnetic field 
strength B = ^B a bdx a A dx b . We also work with the vector densities V a = \e ahc T>ab and 

r>a 1 ,abc R 

-D — 2 6 ^ab- 

We substitute (12.1 ft and ( 12. 2ft into ( 12.31) . we find the 1 + 3 dimensional form of the Maxwell 
equations: 

dV = p, dH-V = j; (2.4) 
dB = , dE + B = . (2.5) 

The 3D exterior derivative is denoted by d and the time derivative by a dot. This is the 
premetric (i.e., metric independent) framework of electrodynamics which summarizes the 
Maxwell equations and their physical interpretation. 
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The properties of the medium under consideration are expressed by means of the consti- 
tutive law. With the assumptions of locality and linearity, we have 

1,1 kl 

Hij = = ^ijmnX mn Fkl • (2-6) 

It is convenient to introduce, besides the constitutive tensor density , the tensor x^ k \ 



since the latter is used conventionally in electrodynamics, see 26]. Its symmetries are 
ytjki _ _yj%ki _ _yAjik^ j e ^ Yias 36 independent components. 

The constitutive law can be put into different forms in order to customize it for different 
applications. If we put it into a 6 x 6 form 




C b a B; 



A ba V h a / \ B b 




(2.7) 



then the 3x3 matrices A, B, C, T> can be related to the 4-dimensional constitutive tensor 
density x 1 by 

A ba = x° a0b , B ba = \ e acd he f X cdef , (2.8) 
C\ = \ e bcd X cd0a , = ~ e acd X° bcd • (2.9) 

Then 6x6 form of x^ kl can be written as 

X IK =( BoiV \ I , (2-10) 
\C\ A ab 

with J, K, ... = 01, 02, 03, 23, 31, 12. 

We can decompose x^ kl irreducibly under the linear group. As we have shown elsewhere 



28), we find 

yAjkl _ {l)^ijkl _|_ (2)^ijkl _|_ (3)^,i jkl _ (2.11) 
36 principal 20 skewon 15 axion 1 

The irreducible pieces carry the additional symmetries 

(l) x m = (l) x klij ^ (2) x ijkl _ _{2)^klij ^ (3) x ijki = (3) x [ijkl] _ (2.12) 

The principal part with its 20 independent components is the only one discussed conven- 
tionally. The skewon part with its 15 components vanishes if one assumes the existence of a 
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Lagrangian 4-form from which the constitutive law can be derived completely. Finally, the 
axion piece with only 1 independent component is totally antisymmetric: 

(3) x y« = ae iiw. (2.13) 

The a is a 4D pseudoscalar. 

We take care of the irreducible decomposition and evaluate the matrix elements of (12.71) : 

tta = {^ab ~ £a bc m c ) B b + {- f3 b a + s a b - 8 b a s c c ) E b - a E a , (2.14) 

V a = ( £ ab _ e abc ^ ^ + ( ^ + _ ga^cj B b + a B a . (2.15) 

We have e ab = e ba , yU^, 1 = fj, ba , and (3 C C = 0. Thus we have the independent components of 
e ab (6), (6), P a b (8), m c (3), n c (3), s a b (9), and a (1). This adds up, as it is required, 
to 36. 

In this paper, we assume that there exists a Lagrangian from which the constitutive law 
can be derived. Thus, ( 2 )^ l J' w = 0. Moreover, we assume a vanishing axion part a = 0. Still, 
in certain substances an axion part can be present, as we discussed recently, see 34, 3s| - 
After these "amputations" , the constitutive law to be investigated reads 

H a = ^B b -f3 h a E h , (2.16) 

V a = e ab E b +(3 a b B b . (2.17) 

Recall that f3 c c = 0. 

III. CONSTITUTIVE RELATION FOR MAGNETOELECTRIC MEDIA 

Throughout the paper we will use the exterior calculus which proves to be very effective 
and convenient both in describing the general formalism and in the specific computations. 
Here we put the constitutive relation for the magnetoelectric media into a simple and trans- 
parent form by using the language of exterior calculus. 

As is worked out in 28(, magnetoelectric properties of the medium are described by 
a tracefree 3x3 matrix (3 a b in the rest frame. The corresponding spacetime foliation is 
called the laboratory foliation, with the coordinate time variable a labeling the slices of this 
foliation. The spacetime metric g introduces a scalar product in the tangent space and 
defines the line element which reads (a, b, ... = 1, 2, 3) 

ds 2 = N 2 da 2 + g ab dx a dx b = N 2 da 2 - ^g ab dx a dx b . (3.1) 



Here N 2 = g(n, n) is the length square of the foliation vector field n, and dx a = dx a —n a da is 
the transversal 3-covector basis, in accordance with the definitions above. The 3-metric ^g a b 
is the positive definite Riemannian metric on the spatial 3- dimensional slices corresponding 
to fixed values of the time a. This metric defines a 3-dimensional Hodge duality operator -. 



Modern discussion o 
can be found in 26, 



the classical electrodynamics, in particular, using the exterior calculus, 
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Following Feigel, we do not consider the effects of gravity. Accordingly, we are in the 
Minkowski spacetime and a convenient choice of the laboratory foliation is a = t and n a = 
(hence dx a = dx a ). The spatial metric is Euclidean, ^g a b = diag(l, 1, 1), and N = c. 

We introduce the 1-form j3 a = f3 a t,dx b . Then the constitutive relation for the magneto- 
electric medium, provided we assume isotropic permittivity and permeability, reads 



V = eeo^E - (3 a Ae a \B, 
H = — *-B-(3 a Ae a \E. 



(3.2) 
(3.3) 



We can decompose the tracefree f3 a into two irreducible parts (symmetric, and antisymmet- 
ric): 

^ = {l)pa + (2) = (1)^ + *0 ^ ^ _ ^ 

The antisymmetric (pseudotrace) part is defined by ^(3 a := ~(/3 A dx a ) with a 1-form (3 := 
~(dx a A (3 a ). Obviously, e a \^ (3 a = 0. Finally, the first irreducible part is trace- and 
pseudotrace-free, i.e., e a \^P a = and dx a A W (3 a = 0. 

The first term in (13.41) describes the tracefree symmetric part of the matrix f3 a t>, whereas 
the second term is its antisymmetric part. When the first irreducible piece vanishes, the 
constitutive relation (13.2I) - (I3.3I) becomes much simpler: 



V 

n 



ee -E 
1 



BAP, 
B + ±(EA0). 



(3.5) 
(3.6) 



In this paper, we will mainly analyse the general case fl3.2l) - fl3.3p . with special attention 
to the medium studied by Feigel. The latter is characterized by the magnetoelectric matrix 



/ 



(3 a t 





/3 2 3 



\ 



P : 2 

7 



3 



(3.7) 



Such magnetoelectric susceptibilities can be induced in an ordinary medium (characterized 
by the permittivity e and permeability jj) under the action of external electric and magnetic 
fields applied along the second and third axes. In spite of the simple form of (13. 7p . there are 
still two nontrivial irreducible parts. When f3 2 ^ = j3 3 2, the pseudotrace $ vanishes, whereas 
for j3 2 3 = —/3 3 2 the first irreducible piece disappears. 

Our subsequent analysis reveals that different irreducible parts of the magnetoelectric 
matrix are differently involved into the possible Feigel effect. 

The magnetoelectric matrix has the dimension [(3 a b\ = [a/^o/Vo]- Accordingly, it is 
convenient to introduce a dimensionless object ft \ := /3 a b/A where A = ^eeo/fif^o- We will 
use the same "overlined" notation also for the various exterior forms constructed from the 
magnetoelectric matrix. 

IV. BIREFRINGENCE IN MAGNETOELECTRIC MEDIA 

The Fresnel approach (geometric optics) to the wave propagation in media and in space- 
time with the general linear constitutive law (12. 6p gives rise to the extended covariant Fresnel 
equation [28j for the wave covector g,: 



Here the fourth order Tamm-Rubilar (TR) tensor density of weight +1 is constructed from 
the constitutive tensor x^ M as a cubic contraction with the Levi-Civita densities: 



gijkl 



(x) QiQjQkQl = . 



(4.1) 



gijkl 



(X): 



1 



,mnr(i ~ / j\ps\k l)qtu 





mnpq 



^rstu ~X_ 
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Using the (1 + 3)-decomposed representation in terms of the 3 x 3-matrices (I2.10p . the 
independent components of the TR-tensor H4.2[) read explicitly as follows: 



M : 
M a : 

M ab. 



6? 0000 = det^, 

4 g 000a = -hcd (A ba A ce C d e + A ab A ec V d ) , 

Q g00ab = }_ A (ab) [(^2 + ^3 _ + (C c + V c 

+ (C d c + V d ){A c{a C b) d + V d {a A b)c ) - C d d A c(a C b) c 
-V c {a A b)c V d d - A dc C {a c V d b) + (A (ab) A dc - A d{a A b)c )B dc , 



(4.3) 
(4.4) 



M 



abc 



4£ 



Oabc c de(c\ 



[B df {A ab) Vj -V e a A b)s 



j^abcd ._ gabcd _ ^5{c^gh\d g 

Then, in (1 + 3)-decomposed form, the extended Fresnel equation (14.11) reads 



+ B fd {A ab) C f e - A fla C b) e ) + C a f P e 6) Vj + V f a C b l C f d 

'1 

2 



■i abed 



■ A ab) B ge - C a e P g 6) 



(4.5) 

(4.6) 
(4.7) 



qlM + g 3 g a M a + q 2 q a q b M ab + q q a q b q c M abc + q a q b q c q d M abcd = 



abc 



rabed 



(4.8) 



In Minkowski spacetime, for the general linear constitutive relation (I3.2p - fl3.3l) we have 
explicitly 

A ab = -ee S ab , B ab = —6 ab , C\ = V b a = (3\. (4.9) 

Wo 

In the special case (I3.5p - (I3.6I) when only the irreducible antisymmetric (pseudotrace) part 
of the magnetoelectric matrix is nontrivial, we analysis of the Fresnel equation (14. 8p reveals 
the clear birefringence effect. Indeed, using the constitutive relation fl3.5H - fl3.6j) in (14. 3ft - 
(14. 7p . we find explicitly (with n = ^Jeji as the usual index of refraction, and c = l/^/^oMo 
the vacuum speed of light): 

2 



M 



n\ 3 



M ab = _ \s ab (2 + /3 2 )-5/3 a /3 b ] , 
c L J 

M abc = -25 {ab p c) (2 + p 2 ) + 2~p a p b p c 



M 



abed 



n 



[-5 {ab 5 cd \l + f3 2 )+5 {ab f3 c f3 d) ] . 



(4.10) 
(4.11) 
(4.12) 
(4.13) 



Here from the magnetoelectric 1-form $ we extract the dimensionless covector with the 
components f3 a := e a J/3/A. Further, we denote [3 a := 5 ab (3 b and [3 2 := [3 a (3 a . 
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Substituting this into the Fresnel equation (14.11) . (14. 8p . we find birefringence: the quartic 
wave surface is factorized into the product of the two light cones, 



Q m {x) mmm = - (<?rW) (<? 2 fc W) = o. 



(4.14) 



The two optical metrics 
to 



311 ] depend explicitly on the magnetoelectric properties according 

(4.15) 



9? 



and 



92 



n 2 
c 2 






-S ab 



(4.16) 



-5 ab (l + [3 2 ) +p a p b 

respectively. It is interesting that the magnetoelectric covector manifests itself as an effective 
"rotation" of the spacetime related to the off-diagonal components of the optical metric. 

In the general linear case of the constitutive relation (I3.2p -( l3~3l . the quartic Fresnel surface 
does not factorize into the product of the two light cones, in general. The birefringence effect 
is much more nontrivial in this case. We, however, will eventually specialize to the medium 
(13. 71) studied by Feigel, and in this case we can analyse the propagation of the planes waves 
to the end. 



V. PLANE WAVES IN A MAGNETOELECTRIC MEDIUM 



A general discussion of the propagation of electromagnetic waves in media with a local 



and linear constitutive relation can be found in [4j, |28|, |32|, l33j . Various special cases of the 
constitutive tensor were analysed, including also the magnetoelectric cases. Here we will 
confine our attention to the specific form of the magnetoelectric matrix (13. 7p . 

In order to clarify the possible Feigel effect, we need to analyse not only the waves inside 
the sample (as was done originally in [6?]) but also the waves in the outside vacuum space. 
The appropriate qualitative picture is as follows (see Fig. [TJ): Let us put the magnetoelectric 
matter between the two parallel planes Si = {x = —£} and S 2 = {x = +£}. Feigel considers 
the case when the magnetoelectric properties are induced by external electric and magnetic 
fields. When these external crossed electric and magnetic fields are applied parallelly to the 
boundaries Si and S 2 , the magnetoelectric matrix will have the form (13. 7p . The idea of 
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y 




FIG. 1: An infinitely extended (in y and z directions) slab of a magnetoelectric medium is located 
between the plane boundaries S\ and S*2. We study the propagation of electromagnetic waves in 
x-direction inside and outside the magnetoelectric medium. 

Feigel is that the virtual electromagnetic waves (vacuum fluctuations of the electromagnetic 
field) propagating in such a magnetoelectric medium could produce a nontrivial momentum 
along the distinguished axis x. However, the presence of similar virtual waves outside of the 
medium was not taken into account (despite the fact that the magnetoelectric sample was 
assumed to have finite size). Here we reanalyse carefully the picture, taking into account all 
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TABLE I: Amplitudes in the different regions and direction of the corresponding wave. The polar- 
ization index 1 and 2 describes a wave with the electric field in y- and in z-direction, respectively 



region 1 


region 3 


region 2 




Pl,P2 


Cl,C 2 


h,b 2 


<?1,<?2 


di,d 2 



the virtual electromagnetic waves inside and outside of the medium alike. 

We begin by noticing that outside of the matter, we have the "bath" of the virtual 
photons some of which will penetrate the interior of the sample, reflecting and refracting 
at its boundaries. Obviously, the largest contribution to the possible effect should come 
from the electromagnetic waves which travel along the e., with the wave vectors 

normal to the boundaries. Clearly, for each right-moving wave falling on the left boundary 
Si, there exists an equal but opposite left-moving wave falling on the right boundary S 2 . 
The contributions of these ingoing waves to the momentum density of the electromagnetic 
field are equal with opposite sign, thus providing a balance of the light pressures in the 
left and in the right vacuum regions. However, we have to find the contributions of the 
outgoing waves. If they turn out to be different in the left and in the right vacuum regions, 
this would seemingly yield a violation of the momentum balance and would encompass a 
nontrivial Feigel effect. 

There are three regions (see Fig. [I] and Tabled]): 1) the left vacuum space (for x < —£), 
2) the right vacuum space (for x > £), and 3) the interior region filled with magnetoelectric 
matter (for — I < x < £). The configurations of the electromagnetic field in these three 
domains read, respectively, as follows, where we use the complex representation to simplify 
the formulas: 
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1) In the first region (x < —I): 

E = e~ iwt [( ai e ikx + b ie ~ ikx ) dy + (a 2 e ikx + b 2 e- lkx ) dz] , 
B — — e~ iujt dx A [( ai e ikx - b ie - ikx ) dy + (a 2 e lkx - b 2 e^ kx ) dz] , 

OJ 

V = e e~ iujt dx A [- (a ie lkx + b x e~ ikx ) dz + (a 2 e tkx + b 2 e~ lkx ) dy] , 



E = 




B = 


h 




OJ 


V = 






k 


H = 


e 




/IqOJ 



iujt 



[{a^^ - b x e ikx ) dz - (a 2 e lkx - b 2 e ikx ) dy] . 



(5.1) 
(5.2) 
(5.3) 
(5.4) 



Here the complex amplitudes a\ and a 2 describe the right-moving waves with the two differ- 
ent polarizations (the subscript i refers to the first polarization with the electric field along 
the y-axis, whereas the subscript 2 denotes a second independent polarization with the elec- 
tric field along the z-axis). The complex amplitudes b\ and b 2 describe the corresponding 
left-moving waves. With k = oj/c one can straightforwardly check that this configuration is 
a solution of the Maxwell equations. 

2) Similarly, in the second region (x > £): 



E = e~ wt [( Cl e lkx + d ie ~ ikx ) dy + (c 2 e ikx + d 2 e~ ikx ) dz] , 
B — — e~ iujt dx A [( Cl e lkx - d x e' lkx ) dy + (c 2 e ikx - d 2 e~ lkx ) dz] , 



e -iu,t ^ c±e ikx _ dl6 -ikx^ dz _ ^ikx _ d2( ,~ikx^ dy ^ _ 



E = 






k 


B = 


— e 




OJ 


V = 


E e~ 




k 




H = 







(5.5) 
(5.6) 
(5.7) 
(5.8) 



Now the amplitudes Ci,c 2 and d±,d 2 describe the right- and the left-moving waves, respec- 
tively, with the two polarizations, and again k — uj/c. 

3) In the third (interior) region (with — i < x < I), the field configurations look somewhat 
more nontrivial: 



E 
B 
V 



= e 



—itot 



pie ikfx + ifcj x 



) dy + (. 



p 2 e ik t x + q2e ik 2 x 



dz 



(5.9) 



-iuit 



dx A 



OJ 



= e 



dx A 



(ktPie iktx - Kq ie - ik i x ^ dy + (ktp 2 e ik * x - k 2 q 2 e- ik * x ^ dzj ,(5.10) 



ee - ^/3 2 3 J qi e~^ x dz + (ee + ^\ J p 2 e~^ x dy] , 

(J?L - p\) pie ^ x dz - ( + P 3 2 ) p 2 e^ x dy 

/ k. 



ee + — i 2 , ) Pl e^ x dz + ( ee - ^ >' :i > ) p 2 t 'W/y 
oj J \ OJ 

An 



H = e 



—iujt 



(5.11) 



- (-^- + P\ ) qie~^ x dz + 
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- f?2 g 2 e 



dy 



(5.12) 



The matter is characterized by the permittivity e, the permeability fi, and the magneto- 
electric matrix (3 a b (the latter is given by (13.71) ). The right- and the left-movers are now 
described by complex amplitudes Pi-,p 2 an d <?i, <?2, respectively. The birefringence, however, 
is manifest in the inequality of kf ^ k 2 : the waves with different polarization have different 
propagation vectors. Explicitly, we find 



c 



(/3 2 3 ) 2 + l±/? 2 3 



fc± - — 

2 c 



(5.13) 
(5.14) 



As before, we use here the dimensionless magnetoelectric variable (3 a b '■= (3 a b/^- 

The twelve amplitude coefficients ai, a 2 , &i, b 2 , ci, c 2 , d±, d 2 ,p 1 ,p 2 , Qi, 92 are n °t arbitrary 
but are related among themselves via the jump conditions for the electromagnetic field 
strength and the excitations at the boundaries Si and So- There are, as usual, twelve jump 



conditions - six for every boundary surface. They read 

Au = 0, 



2Sj: 



fay -v<n) 

(B {1) - B {3) ) 
(B {3) - B {2) ) 



Si 



Si 



A v = 0, 



Av = 0, 



5 2 



A i/ = 0, 



s 2 



t a\ (H(i) - W( 3 )) 

r A j - £ (3) ) 

r ^J (W(3) - W( 2 )) 
t aJ (-B(3) - #(2)) 





= 0, 


Si 




= 


Si 






= 0, 


s 2 




= 


S 2 





(5.15) 
(5.16) 
(5.17) 
(5.18) 



Here v = dx is the 1-form density normal to the surfaces and T\ = d y , r 2 = d z (A = 1, 2) are 
the two vectors tangential to the boundaries. 

Substituting f l5.ip -( l5TT2|) . we find that some of the jump conditions are trivially satisfied 
since B A v = and D A v = in all the three regions. Accordingly, we are left with only 
eight conditions which result from the continuity of ta\TC and ta\E at the two boundaries. 
After some algebra, noting in particular that 



es ± ^/3\ 



,± 



with the abbreviations 

Oil 



V 



ee Q =f — /3 3 2 

fiflQUJ 



V A*o ^ 

V /4) 



(5.19) 
(5.20) 



1 [1 + (/9 2 3 ) 2 ], 



«2 := 



(/^ 3 2) 2 ] 



(5.21) 



14 



we can bring this system into the form of the eight equations: 



e iM ai - e ik % - a x e~ 



+ e ik % - e~ ik1[ 



Pi 



Jk7£ 



Jk~£ 



e iM Ci - e iM d x - a x e ik ^ l px + cxye. ikl l 



e iM ci + e~ iu ' 



e a.2 — e o 2 



e iM c 2 - e M d 2 - a 2 e 1 ^ l p 2 + a 2 e 



ik 2 £ 

e ike c 2 + e~ iM d 2 - e ik ? i p 2 - e~ ik2 1 



e~ iM a 2 + e ik % - e~ ik * € p 2 



ik(. 





= 0, 


(5.22) 


Qi 


= 0, 


(5.23) 


<h 


= 0, 


(5.24) 


HI 


— n 


(5 95") 


<h 


= o, 


(5.26) 




= 0, 


(5.27) 




= 0, 


(5.28) 


(h 


= 0. 


(5.29) 



Thus, we can always choose 2 waves (inside or outside) the medium as independent and find 
all the other waves in all three regions of space from the system (I5.22p - fl5.29p . Obviously, the 
waves with a particular polarization are only related to the waves of the same polarization. 
They do not mix with the modes of a different polarization. 

The following three cases exhaust all possible situations: (i) we choose as primary the 
waves inside the medium, i.e., the amplitudes Pa,Qa (with A = 1, 2) are independent, and 
the amplitudes a a, ^a and ca, ^a outside the matter are obtained from them as secondary, 
(ii) the waves in one vacuum region (for example, in the first one) are primary, then a a, &a are 
independent and the waves in the matter P a, Qa and in the second vacuum region ca, dA are 
derived from them, (iii) the ingoing waves in the vacuum regions, ciA,dA, are independent, 
then the waves in the matter P a, c[a and the outgoing waves 6^, ca are derived. The system 
fl5.22p - fl5.29p can be straightforwardly solved for all these cases. 

(i) Assuming p A , q A to be the independent variables, see Table [Til we fi n d from (I5.22p - 
( 15.29P the amplitudes of the waves in the two vacuum regions: 

JH r 

:i + a A )e- ik ^p A + (1 - «aK 



a A 
b A 
ca 



Jk - e QA 



-ikl 



[l-a A )e-^p A + {l + a A y k A l q A 



ikC. 



Akl 



1 + a A y k+Al p A + (1 - a A )e 



-%k~£. 



dA = —r- 



[l - a A )e ik ^p A + (1 + a A )e- ik ^q A 



(5.30) 
(5.31) 
(5.32) 
(5.33) 
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TABLE II: Case 1. Given are the amplitudes in region 3, the rest is computed 



region 1 


region 3 


region 2 


compute 


Pl,P2 


compute 


compute 




compute 


TABLE III: Case 2. Given are the amplitudes in region 1. The rest is computed. There is also an 
equivalent case if only the amplitudes in region 2 are specified 


region 1 


region 3 


region 2 




compute 


compute 




compute 


compute 



(ii) Assume now that the waves in the left vacuum region are primary, see Table IIII1 
Then a^, &a are independent variables, and we find for the amplitudes in the medium and 



TABLE IV: Case 3: The amplitudes of the ingoing waves (with respect to region 3) are specified, 
the rest is computed 



region 1 


region 3 


region 2 


ai,a 2 


compute 


compute 


compute 


compute 


di,d 2 
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in the second vacuum region: 



PA 



ca 



d,- 



2a A 
2a A 
2a A 



[(I + a A )e- M a A -(I- a A )e M b A ] 



2a A 



[-(l-a A )e- iki a A + (l + a A )e ike b A ], 

[K A e~ 2lkl a A - z(l - a\) sm(k+ + k A )tb A ] , 
[i(l - a 2 ) sin(£;+ + ^)£e" 2lfc V + ■ 



Here we denoted 



JsTa := 2a A cos(k\ + k A )£ + i(l4-a A )sm(k\ + k A )£, 
A A := i^if* = Aa 2 A + (1 - c^) 2 sin 2 (A;+ + k A )t 



(5.34) 
(5.35) 
(5.36) 
(5.37) 

(5.38) 
(5.39) 



The second quantity will be needed below. The star denotes complex conjugation as usual. 

(iii) Finally, if we assume, in accordance with Table HIl that the ingoing waves in the 
two vacuum regions, namely a A ,d A , are independent, the amplitudes of the waves in matter 
p A , q A and of the outgoing waves b A , c A read: 
K A e~ iM 



PA 

b A 
c A 



A A 

K A e~ iM 

K A e~ 2iM 
K A e~ 2iU 



1 + a A )e~ ik ^ e a A - (1 - a A )e ik ^ e d A 
1 + a A )e~ ik A l a A - (1 - a A )e ik ^ e d A 
-i(l - a\) sm(k\ + k A )£a A + 2a A e i ^~ k ^ )e d A 
2a A e l{k ^- k ^ e a A - z(l - a\) sin(£;+ + k A )£d A 



Both quantities (15. 38ft and (15.391) depend on the sums 



kt+k A =—SJVi, 



(5.40) 
(5.41) 
(5.42) 
(5.43) 

(5.44) 



with fix = /3 2 3, f3 2 = — /3 3 2- In practice, the magnetoelectric parameters are rather small 



(typically of order of 10 4 — 10 6 ), so with very high accuracy, we have k A + k A 



2nu) 



2k. 



However, since the difference 



k A k A 



2nu — 

PA 



(5.45) 



for the magnetoelectric matter is nontrivial, the amplitudes of the left-moving waves in 
matter are in general distinct from that of the right-moving waves. Similarly, the amplitudes 
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of the outgoing waves in the two vacuum regions are different, in general. We have to check 
now if such a difference can be manifest in different field momentum densities and forces in 
the two vacuum regions. 



VI. ENERGY, MOMENTUM, AND FORCES 



In vacuum, the energy-momentum 3- form of the electromagnetic field reads 28] 



S Q = ~ [F A (e a \H) - (e a J F) A H] . (6.1) 

Using the (1 + 3)-decomposition, with F = E A dt + B and H = —H A dt + V, we find 
explicitly for the temporal and spatial parts: 

S = u-dtAs, (6.2) 
S a = -p a - dt A S a . (6.3) 

As we see, the energy-momentum generically decomposes into the four pieces: The energy 
density 3-form 

u := ^(EAV + B AH) , (6.4) 
the energy flux density (or Poynting) 2-form 

s:=EAH, (6.5) 

the momentum density 3-form 

p a := - B A (e a \V) , (6.6) 
and the Maxwell stress (or momentum flux density) 2-form of the electromagnetic field 

S a := \ [(e a \E) AV-(e a \V) AE 

+ (e a \H) AB-(e a \B) AH] . (6.7) 

We will evaluate these expressions for the plane wave configurations discussed in the 
previous section. Only the mean averaged (over a time period) quantities have a direct 
physical meaning. Using (I5.ip - fl5.4l) . we then find for the averaged quantities in the first 
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region: 



<u> 



<s> 



<Pa> 



y ^(|a A | 2 + \b A \ 2 ) dxAdyA dz, 
2 ) dy A c?2;, 



/ 

2 



^0 \ ^ 



(|ou| 2 - |6a| 2 )^ Ady Adz 





/ 



(6.8) 
(6.9) 

(6.10) 



<5'a> 



£o 
2 



/ 



V 



(\a A \ 2 + \b A \ 2 ) dy A dz 





(6.11) 



/ 



The form of corresponding quantities in the second region is the same, with the amplitudes 
a a replaced with ca and b A replaced with d A . 

Consequently, it remains to calculate the moduli of the amplitudes of the incoming and 
outgoing waves, and compare the resulting quantities in the first and in the second regions. 
We begin with the case (i) which directly corresponds to the computations of Feigel. We take 
the virtual waves in the medium as primary fields, assuming the equal amplitudes pa = q A 
for the left- and right-movers, and then find from (I5.30p - fl5.33p the waves in the two vacuum 
regions: 

\Pa\ 2 



°A 



\d, 



[1 + a 2 A + (1 - a 2 A ) cos(k\ + k^)£] . (6.12) 



Accordingly, we find that the total momentum and the total energy flux are vanishing in 
both vacuum regions, see (16.91) and (I6.10p . whereas the energy density and the stress densities 
are the same, see f)6.8p and (16.111) . Recalling that the force acting on the boundary can be 
calculated as an integral T x = — J <S X >, we then conclude that the resulting force acting 
on the medium is zero, since on the left boundary the normal vector points in the negative 
x-direction, and on the right boundary in the positive one. 

The same conclusion is derived when we analyse the case (ii) with the waves of the equal 
amplitude travelling in the first or in the second vacuum regions. When we evaluate the 
combined effect of their contributions, we again find the zero resulting force acting on the 
medium. 
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Somewhat different is the case (iii) when the ingoing waves (the right-moving oa in the 
first (left) region and the left-moving wave in the second region with the same amplitude 
d,A = a a) are considered as primary. Then for the outgoing waves we find from (j5.42j) - (l5.43p 

4a A , 



M 2 = 


\ a A\ 


\ca\ 2 = 


M 



1 + 
1 - 



A A 

Aoa 



'1 — a A ) sin(k A + k A )£ sm(k A — k A )£ 
[1 — a A ) sin(k A + k~ A )£ sin(fc^ — k A )£ 



(6.13) 
(6.14) 



As we can see, the contributions of the outgoing waves to the field momentum are clearly 
different in the two vacuum regions. The difference reads explicitly 

8\a A \ 2 a A {l - o? A ) 



v>a\ - M 



A, 



sin(fc^ + k A )£ sin(A;^ — k A )£. 



(6.15) 



When (3 a b = 0, in view of (I5.45P the "bath" of virtual waves around the sample is in 
equilibrium since then (I6.15P vanishes. Under this assumption, the total momentum of the 
waves in both vacuum regions is obviously equal to zero. However, for magnetoelectric 
matter, the mentioned "bath" is still balanced in the sense that the field momentum carried 
by the waves in the left vacuum region is the same as that of the waves in the right vacuum 
region. Namely, substituting (16.131) and (16.141) into (16.101) . we find that the momentum 
densities of the electromagnetic field in both regions are equal 



<Px>- 



2 |^ |2q, 

(1 — a\) sin(A;^ + k~ A )£ sin(A;^ — k~ A )£dx Ady A dz. (6.16) 



A=l 



cA, 



At the same time, the stress is different in the two regions. When we compute the cor- 
responding force T x = — J < S x >, acting on the boundary of the sample, the resulting 
expressions will read for the first (left) and for the second (right) surfaces, respectively: 

2 . . T 9a „ 1 

(6.17) 



J^jf*" = ^ Sq I Oa 1 2 *A 
A=l 

2 

-r-right I 1 2 a 



I + ^(1 - a\) sin(A;+ + k~ A )£ sin(£;+ - k A )£ 



A=l 



A A 

_ 2aA i 
A A 



1 — a A ) sin(fc^ + k A )£ sin(A;^ — k A )£ 



(6.18) 



Here A is the area of the boundary surface (we assume the left and right surfaces to be 
equal). Thus, there will be a nontrivial resulting force acting on the sample in the direction 
of the magnetoelectric vector: 



jrleft + fright _ 



A=l 



A,. 



sm(k A + k A )£ sm(k\ — k A )£. 



(6.19) 
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At first sight, the results obtained, namely (16.161) and (I6.19p . provide a theoretical support 
for the possible Feigel effect. For completeness, however, it is necessary to analyse also 
the situation when the directions of all waves are reversed, i.e., instead of assuming equal 
incoming waves, we should also consider the case of equal outgoing waves. Fortunately, it 
is not necessary to perform a new computation. All we need is to put k — > —k in (15. ip - 
(I5.12p and then note that in the jump equations fl5.22p - fl5.29p we have to change the sign 
of a a — ► —otA- Then repeating the computations of the amplitudes pa, Qa, oa, ca, we arrive 
again to the solution fl5.40p - fl5.43p with the replacements k — > —k and — > —oca- As 
a result, the total momentum density turns out to be again (I6.16p . [It is important to 
note that here we do not have to replace a a — » — «a, since Qa^l 2 — |&a| 2 ) is changed to 
{\°a\ 2 — |°a| 2 ) in (16.101) . and analogously, (\ca\ 2 — \oIa\ 2 ) is changed to {\oIa\ 2 — \ca\ 2 ) hi the 
similar formula in the right region]. However, the resulting force computed from the stress 
on the left and right boundary surfaces will have the opposite sign (and equal magnitude) 
to that of ( 16.19p . Correspondingly, when we consider both contributions together, the total 
force will be found to be equal to zero. In other words, in contrast to Feigel's result, the 
magnetoelectric body will not move, despite the presence of a certain asymmetry between 
the left- and right-moving waves in the matter. 

Several remarks are in order. The above conclusions are based on the evident symmetry 
that characterizes the "bath" of the virtual photons ("vacuum fluctuations") in the regions 
1 and 2: for each left-moving virtual photon there is an equal right-moving virtual photon. 

Furthermore, our observations demonstrate the difference between the two irreducible 
parts of the magnetoelectric matrix /3 a b- Namely, assuming that only first (symmetric) 
irreducible part is nontrivial, we have /3 2 3 = /3 3 2, and then the waves of the different 
polarization contribute to the above formulas with terms of the opposite sign since then 
sin(/c ] f — k^ )£ = — sin(A;^ — k^)L The natural assumption that the virtual waves of any po- 
larization are produced with equal probability then leads to the conclusion that the net result 
will be zero due to the mutual compensation of the contributions of the waves of different 
polarizations. In other words, we find that the first irreducible part of the magnetoelectric 
matrix is irrelevant for the possible Feigel effect. However, for the second (antisymmetric) 
part the situation is different. Then /3 2 3 = — /? 3 2 , and the waves of both polarizations now 
contribute with terms of the same sign. This was also observed in the original computation 
{(J with the result obtained there which is proportional to the skew-symmetric combination 
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(3 2 3 — /3 3 2- This shows that the only the second (antisymmetric or pseudotrace) irreducible 
part of the magnetoelectric matrix (i.e., the 1-form p) is responsible for the possible Feigel 
effect. 

Finally, we can make some predictions going beyond the scheme of the virtual waves of 
the original Feigel picture. Let us assume that real electromagnetic waves are directed on the 
magnetoelectric sample from the two sides. This can be easily achieved if we split the initial 
beam into the two beams which are then, with a simple system of mirrors, directed on the two 
opposite sides of the magnetoelectric sample. Then we are exactly in the situation describe 
as the case (iii) above. The outcome which we derived in this case is that the sample will be 
affected by a nontrivial force (I6.19P from the electromagnetic fields falling on it. The sign and 
the magnitude of this force is determined by the polarization of the infalling waves and by 
the magnetoelectric parameters (3 2 3,(3 3 2 of the medium. This effect can be directly verified 
experimentally. Moreover, one can use this real Feigel effect for the measurement of the 
magnetoelectric parameters (3 2 3,(3 3 2- This method is drastically different and much simpler 
in its practical realization than the usual measurements of the magnetoelectric parameters 
from the observation of magnetization of a medium in an external electric field (resp., electric 
polarization of a medium in an external magnetic field). 



VII. DISCUSSION AND CONCLUSION 



In this paper we analysed the propagation of electromagnetic waves in a magnetoelectric 
medium. Our aim was to use this analysis for the discussion of the possibility of the Feigel 
effect [g] that was predicted recently. 

We derived a natural decomposition of the magnetoelectric matrix /3 a b into the two ir- 
reducible parts: the first (tracefree symmetric), and the second (antisymmetric, or pseu- 
dotrace). The magnetoelectric matrix is tracefree. However, the axion also describes the 
magnetoelectric effect, see (I2.14p and (I2.15p . Being a 4D pseudoscalar, the axion can be 
formally treated as a trace of generalized magnetoelectric 3x3 matrix. We discussed the 



physical properties of the axion piece elsewhere [34j, [35] and here we assume (like in [6| and 
in the previous work 25j) that the axion piece is absent. Ultimately we find that only the 
second part (i.e., the magnetoelectric pseudotrace or antisymmetric part) can be responsi- 



ble for the possible Feigel effect. Physically this result appears to be quite natural because 
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the magnetoelectric pseudotrace 1-form /3 selects a distinguished direction in space. Along 
this direction, we then subsequently find the nontrivial forces and electromagnetic momenta 
acting on the magnetoelectric medium. 

We show that taking into account the finite size of a magnetoelectric sample, one can 
reduce the problem of computation of the electromagnetic energy, the momentum and the 
forces to the vacuum regions just outside the sample. We then demonstrate that the net 
force on the sample turns out to be zero for the careful count of the contributions from all 
the the virtual electromagnetic waves excited in the medium. 

At the same time, we discover that for real electromagnetic waves, that fall along the 
direction of the pseudotrace (3 on a sample from the two opposite sides, is nontrivial. Ac- 
cordingly, it would make the magnetoelectric sample move along /?, indeed, as was predicted 
in the original work (|. We propose to use this observation for the measurement of the 
magnetoelectric susceptibilities. Such an method could provide a useful alternative scheme 
for the measurements of the magnetoelectric parameters, along with the traditional meth- 
ods based on the measurements of polarization induced in an external magnetic field (or of 
magnetization induced by an external electric field). 

..... n 

In the current paper, we have confined our attention (like also Feigel in [a]) to the waves 
induced along f3. Strictly speaking, the complete analysis requires also to take into account 
the waves that move along other directions. The extension of the above results to such 
waves is straightforward, although the formulas (especially those that relate the waves in 
the media and in vacuum at the boundary surfaces) turn out to be appreciable longer. The 
corresponding analysis, however, shows that the main contribution to the possible Feigel 
effect comes from the waves moving along $, so the additional modes do not change the 
conclusions qualitatively. 
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